Low Temperature Behavior of the Vortex Lattice in Unconventional Superconductors 
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We study the effect of tlie superconducting gap nodes on 
the vortex lattice properties of high temperature supercon- 
ductors at very low temperatures. The nonlinear, nonlocal 
and nonanalytic nature of this effect is shown to have mea- 
surable consequences for the vortex lattice geometry and the 
effective penetration depth in the mixed state as measured by 
muon-spin-rotation experiments. 



I. INTRODUCTION 

The presence of nodes in the superconducting gap is 
probably one of the most significant features of high-Tc 
superconductors which has attracted considerable theo- 
retical and experimental attention in recent years. Many 
experiments |l|,g| have confirmed the existence of gap 
nodes. These experiments most commonly indicate an 
order parameter with dx2_y2 symmetry B]. The fourfold 
symmetric nature of the d-wave order parameter together 
with the presence of gap nodes on the Fermi surface open 
possibilities for novel effects to be observable in cuprates. 
An early theoretical investigation of the weak-field re- 
sponse of a dx2_y2 superconductor by Yip and Sauls ||] 
predicted a direction dependent non-linear Meissner ef- 
fect, associated with the quasiclassical shift of the ex- 
citation spectrum due to the superflow created by the 
screening currents. Maeda et al. H reported experimen- 
tal evidence for such an effect in Bi2Sr2CaCu20y, but 
subsequent experiments failed to confirm their findings 
and the situation remains controversial. A similar effect 
was also studied independently by Volovik M. In the 
mixed state, he predicted a contribution to the residual 
density of states (DOS) proportional to the inter-vortex 
distance ^ VH. Such contribution was identified in the 
specific heat measurements on YBa2Cu307^5 (YBCO) 
by Moler et al. g] but later this interpretation was dis- 
puted by Ramirez H who found evidence for similar ef- 
fect in a conventional superconductor VaSi and by others 
[g[. Kosztin and Leggett Q predicted that the nonlocal 
response at very low temperatures will lead to a T^ de- 
pendence of the penetration depth in clean samples in 
contrast to the linear T-dependence obtained from the 
local theory for d-wave materials. 

In the mixed state also, it is conceivable that these ef- 
fects still present themselves in some measurable proper- 
ties such as vortex lattice geometry, magnetic field distri- 
bution etc.. Neutron scattering ||l^ and STM ||ll|] experi- 
ments on the high-Tc compound YBCO revealed a vortex 



lattice structure different from hexagonal - which is ex- 
pected for an isotropic superconductor. However, this 
may be explainable by penetration depth anisotropy and 
twin-boundary pinning without involving any effects as- 
sociated with gap nodes Q . Muon-spin-rotation (^SR) 
experiments |l3-n^, on the other hand, have demon- 
strated an unusual magnetic field dependence in their 
line-shapes for the magnetic field distribution. This has 
been attributed to a field dependent penetration depth - 
which is expected in the Meissner state because of quasi- 
particle generation at gap nodes |^ . It was also modeled 
using an approach based on the Bogoliubov- de Gennes 
(BdG) equations in a square lattice tight-binding model 

At intermediate fields Hd < H <^ Hc2 properties of 
the flux lattice are determined primarily by the superfluid 
response of the condensate, i.e. by the relation between 
the supercurrent j and the superfluid velocity u,. In 
conventional isotropic strongly type-II superconductors 
this relation is to a good approximation that of simple 
proportionality. 
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where ps is a superfluid density. More generally, how- 
ever, this relation can be both nonlocal and nonlinear. 
The concept of nonlocal response dates back to the ideas 
of Pippard iQ and is related to the fact that the cur- 
rent response must be averaged over the finite size of the 
Cooper pair given by the coherence length ^o- In strongly 
type-II materials the magnetic field varies on length scale 
given by the London penetration depth Aq, which is much 
larger than ^o and therefore nonlocality is typically unim- 
portant unless there exist strong anisotropies in the elec- 
tronic band structure |g^. Nonlinear corrections arise 
from the change of quasiparticle population due to su- 
perflow which, to leading order, modifles the excitation 
spectrum by a quasiclassical shift |f9[ 



£k = Ek + Vf ■ Vs, 



(2) 



where Ek = ^e^ -I- A^ is the BCS energy. Again, in 
clean, fully gapped conventional superconductors this ef- 
fect is typically negligible except when the current ap- 
proaches the pair breaking value. In the mixed state 
this happens only in the close vicinity of the vortex cores 
which occupy a small fraction of the total sample vol- 
ume at fields well below Hc2- The situation changes dra- 
matically when the order parameter has nodes, such as 



in dx2_y2 superconductors. Nonlocal corrections to (13) 
become important for the response of electrons with mo- 
menta on the Fermi surface close to the gap nodes even 
for strongly type-II materials. This can be understood 
by realizing that the coherence length, being inversely 
proportional to the gap ||l^, becomes very large close to 
the node and formally diverges at the nodal point. Thus, 
quite generally, there exists a locus of points on the Fermi 
surface where £, :s> Xq and the response becomes highly 
nonlocal. This effect has been first discussed by Kosztin 
and Lcggctt for the Meissner state and by us [p2| in the 
mixed state. Similarly, the nonlinear corrections become 
important in a d-wave superconductor. Eq. (0) indicates 
that finite areas of gapless excitations appear near the 
node for arbitrarily small Vs- This effect has been stud- 
ied in the Meissner state P] but its consequences have 
remained largely unexplored for the mixed state. 

In our previous works pl| , p2[ , we discussed the effect 
of fourfold anisotropics associated with nonlocal response 
on the field and ternperature dependence of vortex lattice 
structure. In Ref. [£l|, starting from a phenomenological 
Ginzburg-Landau (GL) theory, we derived the leading 
fourfold anisotropic corrections to the London equation 
making the usual assumption that the free energy was 
an analytic functional of the order parameter and field. 
We showed that such corrections result in a field-driven 
continuous transition from triangular to square vortex 
lattice. Being derived from a GL theory this London 
equation is expected to lose its validity at low tempera- 
tures and a microscopic theory is needed to address this 
regime. In Ref. |23], we investigated the effect of the 
nonlocality due to the presence of the nodes in the su- 
perconducting gap using a simple weak coupling micro- 
scopic model. At high temperatures a nonlocal correction 
similar to the one suggested in Ref. [^ was obtained. 
At low temperatures however, a novel singular behavior 
was found directly related to the nodal structure of the 
gap which completely changes the form of the London 
equation. This singular behavior has profound implica- 
tions for the structure of the vortex lattice which, as a 
function of decreasing temperature, undergoes a series of 
sharp structural transitions and attains a universal limit 
at T = 0. 

In Ref. |23|, we neglected the effect of the nonlinear 
terms resulting from the excitation of quasiparticles at 
the gap nodes, assuming that they are small compared to 
the nonlocal corrections. In this article, we consider both 
the nonlinear and the nonlocal terms assuming that the 
effects are additive and do not affect each other. We show 
that, as we claimed in Ref. |g^ , the dominant effect which 
determines the vortex lattice geometry and the effective 
penetration depth as defined in /iSR experiments is the 
nonlocal corrections, while the nonlinear corrections play 
a secondary role at low T. 



II. GENERALIZED LONDON EQUATION 

A. Nonlinear Corrections 

Let us for the moment neglect any nonlocality effect 
and focus on the quasiparticles generated at the gap 
nodes. Excitation of quasiparticles at the gap nodes pro- 
duces a current density flowing in the direction opposite 
to the superfluid velocity - sometimes called back-flow. 
The total current can then be written as 



j = -epsVs +jqp 



(3) 



where ps is the superfluid density and Vg is the superfluid 
velocity defined by |Ea 



2 c 



(4) 



A is the vector potential and cf) is the phase of the order 
parameter. The contribution of the quasiparticles gener- 
ated at the nodes to the total current is given by M 



j,p = -4eiV// dsvfis) duwe + ^{sy 

Jfs Jo 



+ Vf{s) -Vs) 



(5) 



where Nf is the density of states at the Fermi surface 
and s is a parameter which represents a point on the 
Fermi surface. A(s) is the superconducting gap which 
in general can have s-wave, d-wave, or other symmetries. 
At zero temperature Eq.(P) leads to 

j,p = -4eNf f dsvf{s)e{-Vf ■ Vs - \A\)^ {vf ■ v^y ~ A^ 



(6) 

At higher temperatures however, this gives the first 
term in the Sommerfeld expansion which is a good ap- 
proximation as long as r < T* = Tc{H/Ho), where Hq 
is of order of the thermodynamic critical field. He [||. 
The presence of the 6*- function in Eq.(|2), results in ex- 
citations only at the nodes in the opposite direction to 
Vs- Fig.n illustrates a circular Fermi surface with a d- 
wave gap. The quasiparticles are excited at the nodes 
marked by Vfi and v/2 in the opposite direction to Vs- 
For a small enough v^, the excitations stay very close to 
the gap nodes. Therefore, one can linearize the gap func- 
tion near the nodes writing A{9) ~ 7A06' with Aq the 
maximum gap and 7 defined by 



^ = K^[^^(^)]-- 



(7) 



The component of jqp along the x'-direction which is di- 
agonal to the X and y (a and b) directions (as illustrated 
in Fig.p is then 




FIG. 1. Circular Fermi surface with a d^2_y2 gap. Quasi- 
particles will be excited at nodes marked by Vfi and w/2 op- 
posite to Vs. 
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eps — cos a cos a 



eps 



Vo 



(8) 



a is the angle between v^ and x'-axis, Ps = Nfv'i is the 
superfluid density, vq = 7A0/W/ is some characteristic 
velocity and 0c is a cut off imposed by the ^-function in 
Eq.(|6|). Similarly the y'-component is 



JqpV' — S-Ps 



Vsy' \Vsy'\ 



Vo 



(9) 



and the total current thereby becomes 



i = -eps[Vs - [Vsx'\Vsx'\x' +Vsy>\Vsy'\y')/vo\. (10) 

The nonanalytic nature of the effect is evident from this 
equation. 

It is now possible to write a free energy in such a way 
that Eq.(|lfl) can be obtained by minimization with re- 
spect to A. Keeping in mind thatj = (c/47r)V x B = 
(c/47r)V X V X A and d/dvg = -d/dA, the correspond- 
ing free energy density can be written as 



f ^ Ps 
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In general, it is possible to solve Eq. ( px| ) for Vg in terms 
of j — (c/47r)V X B, substitute it into Eq.(pl]) and write 
a London free energy only in terms of B and its deriva- 
tives. However, instead of solving Eq.([lO[) exactly, we 
find Vg perturbatively assuming that the nonlinear part 
is much smaller than the linear part. This results in a 
polynomial correction to the London equation which is 



more convenient for the numerical purposes. To first or- 
der in perturbation theory: 



[\J xB- 
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{{W X B),,\W X B\,,x' 
(12) 



Substituting this into Eq. (^l]) and keeping the lowest or- 
der terms, the London free energy density becomes 



h 



1 r 

8^ 



B' + XliX7xBf + A^^mxB) .'3 



^37' Bo 



+ \{VxB)y,\^) 



(13) 



where Aq = y^c^ /Ane^ps is the zeroth order penetra- 
tion depth, ^0 = Vf/nAo is the coherence length. Bo = 
0o/27rAg is a characteristic field of the order of Hd and 
TTc/e is the flux quantum. The first two terms in 



Eq. ( |13D are the ordinary London free energy terms while 
the remaining two terms constitute the leading nonlin- 
ear correction. For magnetic fields in the z-direction, the 
London free energy density becomes 



f 1 



B' + XliVBf + Ak^(\d^,B\' + \dy,B 



37 Bo 



(14) 



and the corresponding London equation will be 
,27r^oA§ 2 



-XoV^B + B-{- 



7 Bo 



{d',,B\d.,B\+d^y,B\dy,B\)=Q 

(15) 



A similar London equation is also derived by Zutic and 
Vails who investigated the effect in the Meissner state 

The most commonly used form for a d-wave gap is 
A{9) = Ao cos(26'). In this case, Eq.(0) leads to a 7 = 2. 
As discussed in more detail below, in order to find the 
magnetic field distribution in a vortex lattice, one has 



to insert a source term Sjp(r 



on the right hand 



side of Eq. (|15| ) with fj being the position of the vortices 
in the lattice. The function p(r) takes into account the 
vanishing of the order parameter at the center of the vor- 
tex cores. Numerically, it is more convenient to work in 
Fourier space rather than in real space. However, because 
of the nonanalyticity of the nonlinear term, the Fourier 
transformation of this term cannot be done by simple 
convolution integrals and numerical techniques such as 
Fast Fourier Transformation (EFT) are required. Fourier 
transforming Eq.(n3) with a proper source term on the 
right hand side yields 



Bk 



Xlk^Bj:-G^i{k,Bi:)^BF{k) 



(16) 



where fc is a reciprocal lattice wave vector, Gni is the 
Fourier transform of the nonlinear term and B is the 



average magnetic field. The cut off function F{k) comes 
from the Fourier transformation of the source term and 
removes the divergences by cutting off the momentum 
sums. 



B. Nonlocal Corrections 

One can add the nonlocal effect to Eq. (|6|) by replacing 
the second term with a nonlocal one as discussed in Ref. 
|P2| . The resulting London equation will be 



B^ + C,jik)hk,Bj: - G„i(fc, Bj:) = BF{k). 



(17) 



Sums over i and j are implicit here. Cij{k) is related to 
the electromagnetic response tensor Q{k) defined in Ref. 
PI by 



/:,y(fc) = Q„(fc)/detg(fc). 



(18) 



The electromagnetic response tensor describes the lin- 
ear response of a superconductor to an applied magnetic 
field. It has a particularly simple form in the gauge where 
V (j) — and w, is proportional to the vector potential A. 
At T = one then obtains ill, 



m 



1 



dsv 



f^Vf3 



2arcsinh y 



yv^+y^ 



(19) 



where y = Vg ■ fc/2A(s). The expression for Q{k) be- 
comes more complicated in an arbitrary gauge and was 
discussed by Millis |g5[. In the present context, how- 
ever, it is easy to verify that, because of the particular 
way Q{k) enters the London equation (|l7|), this is gauge 
invariant as required for an equation containing an ob- 
servable quantity. 

In the local case we have Cij{k) = X^Sij which leads 
back to Eq.([lq). Taking Gni to the right hand side of 
Eq. (p7[) , Bj;: can be obtained by 



% = 



BF{k)+Gni{k,Bf:) 



(20) 



We use this equation to find B^ iteratively for a specific 
lattice geometry. Having Bj\, the free energy can be easily 
calculated using 



^ = -^nl + ^[1 + C,j{k)hk,]Bl 



(21) 



where !F^\ is the free energy due to the nonlinear parts 
as in Eq.(|4l). 



C. Vortex Source 

The functional form of p{r) and thereby F{k) depends 
on the detailed form of the order parameter at the center 



of the vortex and therefore requires a more fundamental 
theory to be evaluated. In the Ginzburg-Landau (GL) 
limit, there have been some conjectures about the pro- 
file of the order parameter at the vortex center and the 
form of the source term Eq-M. Brandt Efl] suggests a 
Gaussian from for the source term 



p{r) = (0o/2^eo)e 



-r'/^il 



which leads to 



F{k) = 



_ -ek'/2 



(22) 



(23) 



Clem P7| on the other hand, assumes the order param- 
eter to vanish as r/^r"^ + ^^ at the center of the vortex 
core near Tc, where ^„ is a variational parameter propor- 
tional to ^0- It is not difficult to show that this form of 
the order parameter leads to a squared Lorentzian source 
term [B2|. 



p{r) 



^0 a 



F{k) = uKi{u) 



(24) 



where Ki is the modified Bessel function. In the extreme 
type-II case, ^t, = v^^o and u — V2^o^- This form of the 
source term shows quite a good agreement with the ex- 
act solution of the GL equation in the vortex state [Q . 
At low temperatures however, the GL equation is not 
applicable. The most commonly assumed form for the 
order parameter near the vortex core which gives a good 
fit to the numerical solutions of the BdG equation pl| ] 
is i/'(^) oc tanh(r/^). The cut off function resulting from 
this form of the order parameter fits quite well with the 
Gaussian form of Eq. ( p^ ) up to high Fourier modes [g2| . 
It is worth pointing out that none of the models for the 
source terms discussed above, consider the d-wave nature 
of the order parameter or other anisotropics which might 
be important for the vortex lattice properties. Recent 
numerical computations within the self-consistent BdG 
theory |Q for a d-wave vortex indicate an order param- 
eter with a relatively weak four-fold anisotropy and an 
amplitude relaxing to its bulk value as ~ 1/r^ far from 
the core even at T = 0. Therefore, a more careful consid- 
eration is required to formulate a reasonable model for 
the vortex core in the London theory. We leave a detailed 
discussion of the vortex core to a separate paper » 
from now on assume the Gaussian form of Eq.(^ 
our vortex source. 



III. NUMERICAL CALCULATIONS 

The numerical calculations are performed by employ- 
ing Eq. (g^) to calculate Bg iteratively. At each step Gni is 
calculated numerically using EFT. In our numerical cal- 
culation, we use Ao = 1400A, k — Aq/^o = 68 and 7 = 2. 
Unlike the Ginzburg-Landau free energy, the London free 
energy can not completely determine the vortex lattice 
by a simple minimization. Instead, one has to impose a 
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FIG. 2. Apex angle /? as a function of magnetic field B at 
T = 0. Circles represent the result of the calculation using 
only the nonlocal corrections |E2[. Squares correspond to the 
calculations considering both nonlinear and nonlocal correc- 
tions. 
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FIG. 3. The effective penetration depth as a function of the 
magnetic field. Circles represent the result of the calculations 
considering only the nonlocal effects whereas squares are the 
result of combined calculation considering both nonlocal and 
nonlinear effects. 



set of source terms located at the position of the vortices. 
The functional form of the source terms does not come 
from the London theory and requires more fundamen- 
tal treatment. The free energy must then be minimized 
with respect to the positions of the vortices (determined 
by the source terms). In general, a 2D lattice can be de- 
termined by four parameters. However, a centered rect- 
angular lattice is energetically more favorable than an 
oblique lattice. On the other hand, the vortex lattice 
spacing is fixed by the average magnetic field B {B k, H 
away from Hd). Thus we are left with two variational 
parameters i.e. the lattice orientation with respect to a 
and b directions and the apex angle [3 - the angle between 
the two basic vectors of the lattice. We therefore find the 
vortex lattice geometry by minimizing T in Eq. ( pl| ) with 
respect to the apex angle [3 for different orientations of 
the lattice. At T = 0, the stable orientation of the lattice 
is aligned with the a and b axes. Fig.H shows the results 
of our numerical calculations for /3 as a function of mag- 
netic field. The upper curve (squares) is the result of 
the combined calculation considering both nonlocal and 
nonlinear corrections i.e. using Eq.(po[) and Eq.(pl|). The 
lower curve (circles) on the other hand, corresponds to 
taking into account only the nonlocal corrections, i.e. ne- 
glecting ^ni in Eq.(pi|) and Gni in Eq.(|20|). As it is clear 
from Fig.pl the difference between the two cases is small 
and about one or two degrees. Therefore, the phase di- 
agram given in Ref. |22] retains its validity qualitatively 
even after adding the nonlinear correction to the London 
free energy. 

We also calculate the effective penetration depth Aos 
for different magnetic fields in almost the same way as it 
is calculated from /xSR data ||l^. In these experiments, 
the /iSR precession signal obtained from the experiment 



is fit to a signal obtained by Fourier transforming the 
magnetic field distribution calculated on a hexagonal vor- 
tex lattice with the same average magnetic field using the 
ordinary London model with some cutoff function. The 
A which provides the best fit to the data is considered 
as Aoff. Here, we calculate Aeff in a different way (but 
similar in spirit) using the fact that in the ordinary Lon- 
don model, for a hexagonal lattice and for a large enough 
field, 
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(B-By^AB^^B^J^Ji 



fc#0 



(H-A2fc2)2 



x-^ib^y: 






(25) 



fc#0 



Therefore, associating all the field dependence of AB^ in 
our calculation with the field dependence of an effective 
penetration depth Aoff, we can define Aeff by 



Aeff 



AB^, 



\AB^ 



(26) 



where ABq is the mean squared value of the magnetic 
field Bo(f) — B obtained by applying the ordinary London 
model on a hexagonal lattice with the same average field 
B and with the penetration depth Aq. 

Fig.0 shows the result of our numerical calculation for 
Aeff using Eq.(p6|). The lower curve corresponds to the 
calculations including only the nonlocal correction. The 
upper curve corresponds to the result of the calculations 
using both nonlinear and nonlocal terms. The effect of 



the nonlinear term to the field dependence of Acfr is al- 
most nothing but an overall shift. Fig.^ exhibits mag- 
netic field distribution n{B) defined by 



n{B') 



A 



(frS[B' - B{f)] 



(27) 



at the average magnetic field B = 5.9T, with A being 
the area of a unit cell. The solid line in Fig.^ repre- 
sents the magnetic field distribution calculated from the 
nonlinear-nonlocal London equation. This line-shape is 
then compared with another line-shape (dashed line) ob- 
tained from an ordinary London calculation but with a 
larger value of Aq. /xSR experiments also produce the 
same kind of line-shape but with some additional broad- 
ening due to lattice disorder, interaction of muons with 
nuclear dipolar fields and finite lifetime of muons. The 
resolution of the magnetic field as a result of this broad- 
ening is SB ^ 1Q~^T. The difference between the solid 
line and dashed line in Fig.^ as well as the double peek 
feature of the solid line is therefore not observable by /zSR 
experiments because of these broadening effects. Thus as 
far as these line-shapes are concerned, it is difficult to dis- 
tinguish a nonlinear-nonlocal effect from a simple shift in 
the magnetic penetration depth in the ordinary London 
model. Fig.^ compares the effect of including both non- 
linear and nonlocal corrections to the London equation 
with the effect of including only the nonlocal term. Com- 
paring the two line-shapes, it is apparent that the effect 
of the nonlinear term is small compared to the nonlocal 
term as was emphasized before. 



IV. DISCUSSION 

The ordinary London equation is not adequate to de- 
scribe all the different properties of a vortex lattice in 
high Tc superconductors especially the properties result- 
ing from the presence of the superconducting gap nodes 
or other anisotropics on the Fermi surface. However, as 
we discussed in this paper and also in our previous publi- 
cations [^ 22 1, a generalized London model with appro- 
priate higher order corrections which take into account 
these anisotropic effects, can still provide a fairly simple 
way to calculate different properties of a vortex lattice. 

In Ref . |gl[ , we established a generalized London equa- 
tion which could describe the structure of the vortex lat- 
tice below Tc down to intermediate temperatures. Our 
results in Ref. pi| were in agreement with the Ginzburg- 
Landau calculations l34|. At low temperatures however, 
the suggested generalized London model ceases to be 
valid because of the nonanalyticities resulting from the 
presence of the nodes in the superconducting gap. In the 
present paper and also in Ref. p2l, we generalized our 
London approach to describe the nonanalytic and also 
nonlinear and nonlocal nature of a d-wave vortex lattice 
at low temperatures. Our numerical calculations indi- 
cate that at both high temperature and low temperature 
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FIG. 4. Solid line: Magnetic field distribution obtained 
from the nonlinear-nonlocal London equation with B = 5.9r 
and Ao = 1400 A. Dot-dashed line: Field distribution ob- 
tained from an ordinary London equation on a hexagonal lat- 
tice with the same B and Aq. Dashed line: Result of the same 
ordinary London calculation but with Ao — 1850 A. 



regimes, the nonlocal parts of the free energy plays the 
dominant role in determining the vortex lattice proper- 
ties. 

The equilibrium vortex lattice geometry exhibits novel 
field and temperature dependence owing to the fourfold 
anisotropic effects expressed by the corrections to the 
London equation (Fig.||, also see Ref. |£l|,|2[). The nu- 
merical calculation of the lattice geometry is rather insen- 
sitive to the details of the vortex cores |pq] . The reason 
is that the details of the magnetic field inside the vortex 
cores mainly affect the self-energies of the vortices. In 
the magnetic fields far below ifc2, the vortex lattice ge- 
ometry is mostly determined by the magnetic interaction 
energy between vortices which is insensitive to the precise 
shape of the core. Therefore, our replacement of the vor- 
tex core by a Gaussian source term should be adequate 
for the vortex lattice structure calculations. 

The effective penetration depth Aos also exhibits field 
dependence at low temperatures as it is illustrated in 
Fig.0. Some important points need to be emphasized 
here: 

(i) The field dependence of Aeff is not linear. The varia- 
tion of Aeff with the magnetic field is faster at lower fields. 
/iSR data is only available for a limited range of the mag- 
netic field and the uncertainty of the experimental data 
is noticeable. Thus, it is hard to judge from the /iSR 
data about the linearity of the field dependence although 
a negative curvature comparable to our result in FigJS is 
evident from the /iSR data for a detwinned YBCO sin- 
gle crystal in Ref. Q . At low magnetic fields where the 
/iSR data is available, the relative variation of the effec- 
tive penetration depth in our calculation is about 7% for 
an increase of IT in the magnetic field which is close to 



0.08 




0.00 



-0.001 



0.000 



0.001 



0.002 



AB/B 



FIG. 5. Solid line: Magnetic field distribution obtained 
from the nonlinear-nonlocal London equation. Dashed line: 
Magnetic field distribution obtained from a London equation 
including only the additional nonlocal term using the same 
parameters. 



7.3% variation obtained from /iSR data for an optimally 
doped and 9.5% variation for a detwinned underdoped 
YBCO single crystal using the same cutoff function as 
Eq.dl) for fitting calculations ||l|,|l6|. 

(ii) More importantly, this field dependence of XeS has 
a predominantly nonlocal origin rather than a nonlinear 
one, contrary to what is generally believed. The contri- 
bution of the nonlinear term to the total (minimized) free 
energy is almost one order of magnitude smaller than the 
nonlocal term. What is more important however, is the 
field dependence and (3 dependence of these terms not 
their orders of magnitude at fixed /3 and B. As we men- 
tioned earlier, we consider /3 as a variational parameter 
which has to be fixed by minimizing the London free en- 
ergy. As can be inferred from Fig.g and Fig.|[ the field 
dependence and P dependence of the nonlinear term in 
the free energy is also smaller than the nonlocal term. 
It is worth noting that in the Meissner state, a linear 
nonlocal term can never produce field dependence in the 
penetration depth (as it is usually defined in that state) 
and therefore a nonlinear term is necessary for such an ef- 
fect 0]. In the vortex state however, a nonlocal term can 
result in a field dependent effective penetration depth. 
To understand this, let's neglect the nonlinear term and 
assume a linear but nonlocal London equation. In that 
case the total magnetic field will be the superposition of 
the fields around individual vortices. The magnetic field 
around an isolated vortex is given by 



B{r) = $0 



d^k 



F{k)e 



ik.r 



{2^TYl+C^J{k)k^k, 



(28) 



where $o is the flux quantum and F{k) is the cutoff func- 
tion resulting from the source term. Cij (k) is defined in 
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FIG. 6. Magnetic field as a function of the distance from 
the center of the vortex (a-direction) for an isolated vortex. 
The solid line corresponds to the nonlocal London equation 
whereas the dashed line represents an ordinary London cal- 
culation with the same value of Aq. Bo = B{r — 0) for the 
ordinary London vortex. 



Eq. (p^ . For small values of fc, Cij{k) ~ ^l^ij. Since 
small k corresponds to large r, one expects an isotropic 
field, similar to the local London case, far away from 
the vortex core. For large values of k however, Cij{k) 
has strong fc-dependence with four-fold anisotropy. Thus 
the closer to the vortex core, the more deviation form 
an isotropic ordinary London single vortex is expected, 
as is clearly shown in Fig.0. At low magnetic fields, the 
vortices are far apart and their magnetic fields overlap 
in regions far away from their cores. The properties of 
the vortex lattice should then be similar to the ordinary 
london hexagonal lattice. As the magnetic field is in- 
creased, the vortices come closer to each other. Although 
the profile of the magnetic field around each vortex re- 
mains unchanged (the nonlinear term is neglected), the 
overlap regions will be closer to the vortex cores and will 
be more affected by the nonlocal term. Therefore, it is 
conceivable that at large magnetic fields, the vortex lat- 
tice properties such as the magnetic field distribution will 
be affected by the nonlocal term in the London equation. 
The magnetic field near the vortex core is reduces by 
the nonlocal term as is clear in Fig.g. This is because 
Cij{k) — XqSij is a positive definite tensor for all k and 
therefore the denominator of the integrand in Eq. (|2^) is 
always larger than the corresponding ordinary London 
one. As a result, the magnitude of AB"^ is smaller for 
the nonlocal case and therefore Acs defined by Eq.([2^) 
tends to be larger. This explains why Acs/Ao is always 
greater than one in Fig.y. Fig.^ exhibits the resemblance 
between a change in the magnetic field distribution due 
to the nonlocal term and due to a shift in the ordinary 
London penetration depth. The slight difference between 
the solid and dashed lines in Fig.Q would be unobservable 



in /iSR experiments as a result of the broadening effects. 
Since no field dependence due to the nonlocal term is 
expected in the Meissner state, the Aoff as defined here 
and also in /iSR experiments is expected to be conceptu- 
ally different from what is usually defined as penetration 
depth in the Meissner state, although they are closely 
related. 

(iii) The magnitude and field dependence of Acff is not 
so sensitive to the apex angle (3. In other words, a few 
degrees change in the variational parameter /3 does not 
modify the magnetic field distribution as much as a vari- 
ation in the average magnetic field does. 

(iv) Calculation of Aoff is rather sensitive to the form 
of the vortex source term. The importance of the source 
term in the calculation of AB^ has already been empha- 
sized by Yaouanc et al. J29| . In Ref. |l3| , Aoff is obtained 
by fitting to the /iSR data using both a Gaussian cutoff 
(Eq. ( p3[)) and also the cutoff function proposed by Hao 
et al. |28y29[| . The difference between the two cases is 
significant and about 30% for the magnitude of Aoff and 
even more (for detwinned sample) for the relative varia- 
tion with respect to the magnetic field. This can explain 
the importance of the source term in calculations of the 
effective penetration depth. 

/zSR experiments ||l^,0 on NbSe2, which is believed 
to be a conventional superconductor, also show a field 
dependence in the effective penetration depth although 
it is much weaker than what is observed in high T^ com- 
pounds. Since there is no node in the superconducting 
gap of these materials, the theory presented in this arti- 
cle can not explain this field dependence. However, since 
the size of the vortex core in these materials is large and 
comparable to the vortex lattice spacing for the magnetic 
fields of experimental interest, it is conceivable that a 
significant effect can come from the cores as it is pointed 
out in Ref. [g9|. Thus, a more careful consideration of the 
vortex core might be necessary in order to have a better 
quantitative explanation of the experimental results [p2[ . 



V. CONCLUSION 

We investigated the effect of the superconducting gap 
nodes on the vortex lattice properties at very low tem- 
peratures by a generalized London approach with higher 
order corrections to the free energy. We found that non- 
local effects, arising from the diverging coherence length 
near the gap node, are predominantly responsible for the 
unusual behavior of the vortex lattice geometry and the 
effective penetration depth. The nonlinear effects asso- 
ciated with the shift of the quasiparticle excitation spec- 
trum play only a secondary role, resulting in small (but 
not negligible) corrections at low T. Contrary to the com- 
mon belief, the effective penetration depth, as defined in 
a /xSR experiment, is not a linear function of the mag- 
netic field and is mainly affected by the nonlocal effects. 
This is in a marked contrast to the nonlinear Meissner 



effect in a d-wave superconductor where the correc- 
tion to Aeff arises strictly from the nonlinear term in the 
London free energy and is linear in the field. 
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